Spin diffusion and the anisotropic spin- 1/2 Heisenberg chain 
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Measurements of the spin-lattice relaxation rate 1/Ti by nuclear magnetic resonance for the 
one-dimensional Heisenberg antiferromagnet Sr2Cu03 have provided evidence for a diffusion-like 
contribution at finite temperature and small wave- vector. By analyzing real-time data for the auto- 
and nearest-neighbor spin-spin correlation functions obtained by the density-matrix renormalization 
group I show that such a contribution indeed exists for temperatures T > J, where J is the coupling 
constant, but that it becomes exponentially suppressed for T <C J. I present evidence that the 
frequency-dependence of 1/Ti in the Heisenberg case is smoothly connected to that in the free 
fermion case where the exponential suppression of the diffusion-like contribution is easily understood. 

PACS numbers: 75.10.Jm, 75.40.Gb 
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I. INTRODUCTION 

Knowledge about the dynamical properties of the spin- 
1/2 Heisenberg chain H = J^- SjSj + i, where J > 
is an antiferromagnetic coupling constant, is the link 
between theory and many experiments on compounds 
which are believed to be good realizations of this model. 
Whereas the static properties of the one-dimensional 
Heisenberg model are well understood based on effec- 
tive low-energies theories and its Bethe Ansatz integra- 
bility, many open problems concerning the dynamical 
properties have only been addressed very recently. An 
important question, much interest has focused on, is 
how the integrability of the pure model as well as small 
integrability-breaking perturbations in any real material 
influence the spin (electrical) 1_5 and heat conductivity. 6,7 
A related question important for nuclear magnetic reso- 
nance (NMR) , neutron scattering and Coulomb drag be- 
tween quantum wires refers to the dynamic spin structure 
factor. A detailed analysis of its lineshape at zero tem- 
perature T and small wave-vector q has been presented 
recently in Refs. 8-10. 

Experimentally, most efforts have been concentrated 
on the compound Si^CuOa which is believed to be an 
almost ideal realization of a one-dimensional spin-1/2 
Heisenberg model with a large in-chain antiferromagnetic 
coupling constant J ~ 2000 K and very small inter- 
chain couplings leading to a Neel temperature Tn ~ 5 
K ~ 0.003 J. Its Heisenberg character is supported 
by measurements of the uniform susceptibility at low 
temperatures which are compatible with a logarithmic 
decrease expected due to marginally irrelevant Umk- 
lapp scattering. n ~ 13 Measurements of the thermal con- 
ductivity have revealed strong spatial anisotropies and 
large parts of the heat current along the chain direc- 
tion have been attributed to magnetic excitations. 14 In 
a pure Heisenberg model the heat current is a con- 
served quantity 15 leading to an infinite thermal conduc- 
tivity Kth- 6 The effect of small integrability-breaking per- 
turbations has been investigated in Ref. 7 and it has 
been found that Kth can remain anomalously large un- 
der certain circumstances. A possible explanation of the 



Sr2Cu03 conductivity data has been proposed in Ref. 16 
arguing that phonon and impurity mediated relaxation 
processes dominate. Another important test of the dy- 
namical properties of this system at small frequencies u 
has been provided by NMR measurements of the spin- 
lattice relaxation rate 1/Ti. Particularly appealing is 
the possibility to separate the contributions from wave- 
vectors q ~ and q ~ n, which are the dominant ones 
for low temperatures, by measuring 1/Ti at inequivalent 
lattice sites with different form factors. 13 ' 17 Theoretical 
studies of the spin-lattice relaxation rate have so far been 
based on the calculation of the dynamical structure factor 
in the framework of low-energy effective theories 18,19 or 
the numerical calculation of imaginary-time correlation 
functions. 20 ~ 23 Recent progress in the calculation of real- 
time correlation functions by the density-matrix renor- 
malization group (DMRG) both at zero temperature 24 
and finite temperature 25 has opened a new and so far 
unexplored avenue to tackle this problem. 

In this article I will focus on the 17 NMR measure- 
ments in Sr2Cu03 13 where evidence for a 9 ~ mode 
with diffusion-like character at finite temperature has 
been found. To test whether or not the spin-lattice re- 
laxation rate behaves indeed qualitatively different in the 
Heisenberg than in the free fermion case, I will consider 
the XXZ-mode\ which interpolates between these two 
cases. In Sec. II the basic theoretical framework to study 
spin-lattice relaxation will be layed out and predictions 
by the Luttinger model discussed. In Sec. Ill the free 
fermion case is considered in detail. The interacting case 
is then analyzed in Sec. IV based on real-time data for 
spin-spin correlation functions obtained by the DMRG 
method applied to transfer matrices. 25 In the last sec- 
tion I discuss and summarize my main conclusions. 



II. BASIC THEORETICAL FRAMEWORK 



The Hamiltonian of the XXZ-cha.m is given by 

N 

H = jJ2[S'jS^ +1 + S]S] +1 + AS^S^ +1 -hS^] , (1) 



where J > is an antiferromagnetic coupling constant 
and h the magnetic field. A parameterizes an exchange 
anisotropy and the model is critical for — 1 < A < 1. By 
Jordan- Wigncr transformation the model can be repre- 
sented in terms of fermionic operators tpj 
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The fermions become non-interacting for A = 0. We 
assume for simplicity that the hypcrfine interaction 



Hhf = 2_^ A r IoS r 



(3) 



experiments, however, we often have the situation that 
T ^> | ft, | so that the effect of the magnetic field on the 
electron spins can be ignored. 

In Sr2CuC>3 , measurements of the spin-lattice relax- 
ation rate have been performed on the copper sites, 17 
which carry the S = 1/2 electron spin, as well as on 
the two inequivalent oxygen sites 0(1) and 0(2) 13 (see 
Fig. 1). The hyperfine interaction drops down rapidly 




0(2) 



FIG. 1: Schematic sketch of the spin chain Sr2Cu03 . 



between the nuclear spin Iq and the surrounding electron 
spins S r , where r is the distance in units of the spacing 
between the electron spins, is isotropic. If the hypcrfine 
interaction is the dominant relaxation process, the spin- 
lattice relaxation rate 1/Ti can be obtained by treating 
Hhf as a perturbation inducing transitions between the 
nuclear levels leading to 26 
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Here uin is the nuclear magnetic resonance frequency 
with u>n <C T in all NMR experiments. By Fourier trans- 
form we obtain 



i-i/|w.)^-ta«») 



(5) 



where the transverse dynamic spin structure factor is de- 
fined by 
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and 



A(q)=J2 ctQrA r 



(7) 



If the hypcrfine interaction (3) is anisotropic, we have 
to replace \A(q)\ 2 - (\A x (q)\ 2 + \A y (q)\ 2 )/2 in Eq. (5). 
In spin-chain compounds there is usually no exchange 
anisotropy, i.e., A = 1. In this case we can replace the 
transverse by the longitudinal dynamic structure factor 
leading to 
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Note, that the spin-spin correlation functions here in 
principle have to be evaluated for finite magnetic field. In 



with distance. It is therefore reasonable to assume that 
only the on-site hyperfine interaction Aq and nearest- 
neighbor hyperfine interaction A\ arc important lead- 
ing to Acu(q) = Jdqcxp(iqr)[A 8(r) + A^r ± 1)] = 
Aq + 2AiCos(q). For the 0(2)-site, on the other hand, 
one finds 



\A 0(2) (q)\ 2 = \A\ 2 
and for the 0(l)-site 
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(10) 



where A, B are material dependent constants. The spin- 
lattice relaxation rate at the 0(2)-site can therefore be 
expressed as 

1/1? = |4 2 Sg» (11) 

and the one at the 0(l)-site as 

im = \\B\*[S?{u;) + Sl'{u)]. (12) 

Similarly, the spin-lattice relaxation rate for the copper 
site is given by a sum of S§ z (u)), Sf z (uj) and S2 Z {oj) but 
with prcfactors which depend on the ratio Ai/Aq. 21 I will 
not consider this case here. 

The low-energy excitations of the Hamiltonian (2) have 
either momentum q ~ or momentum q ~ 2kp 1 with 
Fermi momentum kp = ir/2 in the half-filled case con- 
sidered here. By linearizing the dispersion around the 
two Fermi points and expressing the fermionic operators 
in terms of bosonic ones, a technique termed bosoniza- 
tion, the XXZ-mode\ becomes equivalent (up to irrele- 
vant operators) to the Luttinger model 27 



Hll = - I dx 
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(13) 



Here, v is the spin-wave velocity, <j>(x) a bosonic field and 
n(ai) its conjugated momentum satisfying [<fr(x), Ii.(y)] = 
iS(x — y). For this free boson model the dynamic struc- 
ture can be easily calculated and consists of a uniform 



(q ~ 0) and a staggered (q ~ 7r) part where most of 
the spectral weight is concentrated. 18 At the isotropic 
point, A = 1, the staggered part has been shown to lead 
to 1/Tf ~ ln 1/2 (T /T) at low temperatures where T 
is a scale. 19 The logarithmic temperature dependence is 
a consequence of marginally irrelevant Umklapp scatter- 
ing. The staggered component should completely domi- 
nate the spin-lattice relaxation rates for the copper- and 
the 0(2)-sites and the obtained data 13,17 ' 22 indeed show 
reasonable agreement with this theoretical prediction. 

For the 0(l)-sitc, the form factor (10) leads to a strong 
suppression of contributions from q ~ 7r and contribu- 
tions from q <~ should dominate. The uniform part 
of the dynamic structure factor for the Luttinger model 
(13) at A — 1 is given by 



suggested by Thurbcr et al.. 13 

In the rest of this paper we want to analyze whether or 
not the spin-lattice relaxation rate at the 0(l)-site can 
show such singular frequency dependence under the as- 
sumptions that (a) Sr2CuC>3 is well described by a pure 
Heisenberg model with all other degrees of freedom ne- 
glected, and (b) that the hyperfine interaction causes the 
only important relaxation process. Following Eq. (12) the 
longitudinal auto- and nearest-neighbor correlation func- 
tions will be studied. Although these correlation func- 
tions are not directly related to the relaxation rate for 
A 7^ 0, it is helpful to consider this more general case 
because it interpolates between the exactly solvable free 
fcrmion and the isotropic Heisenberg model we are inter- 
ested in. 



S?{q,u) 



2(l_ e -"/ T ) 
T» w T 



5(w-v\q\) (14) 



2lo 



\q\S(0J-v\q\) 



Here v = Jir/2 is the spin- wave velocity. The spin-lattice 
relaxation rate for the uniform part is then 
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The (5-function peak in the dynamical structure fac- 
tor (14) even at finite temperature is a consequence of 
Lorentz invariance: A single boson with momentum \q\ 
always carries energy u> — v\q\. This simple result will be 
modified by irrelevant operators neglected in (13) corre- 
sponding to band curvature terms. The effect of these 
terms at zero temperature has been analyzed in Refs. 8— 
10. It is, however, not obvious how to generalize these 
results to temperatures T 3> u>. 

The staggered part of the dynamical structure factor 
at A = 1 is given by 18 

nrfi- ;<»-v<i) r ( I - [!£±ai\ 
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where the amplitude D = (27r)~ 3 / 2 has been determined 
in Ref. 28. I ignore multiplicative logarithmic corrections 
here because they are not important for the temperature 
range of interest. Using again the form factor (10) we find 
(1/Tf) s w 0.063 \B\ 2 T 2 for u -* 0. It has been pointed 
out in Ref. 13 that the value for \jT\T predicted by (15) 
agrees with the extrapolation T — > of the experimental 
data but that the slope of 1/TfT in experiment is an 
order of magnitude larger than the one obtained from the 
q ~ 7r contribution in field theory. Most important, there 
is no singular frequency dependence within the Luttinger 
model which contradicts the behavior of the spin-lattice 
relaxation rate at the 0(l)-site 
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III. FREE SPINLESS FERMIONS 

For A = the Jordan- Wigner transformation yields 

S z r z (t)=(S z (r,t)S z (0,0)) (18) 

= ([ C t (r, t)c(r, t) - 1/2] [ C t (0, 0)c(0, 0) - 1/2]) 
= (c+(r,*)c(0,0))<c(r,t)ct(0,0)) 
1 
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where e(k) = — Jcosfc and n k = (c~ J ^ cosfe + 1) _1 . Set- 
ting J = 1, the Fourier transform of the autocorrelation 
function is given by 
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and for the nearest-neighbor correlation function by 
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(20) 



We are interested in the case (3 — 1/T ^> 1, u -C 1 with 
lu <C T. Then the most important contributions to the 
integrals in (19,20) come from e ~ (Fermi points) and 
from e ~ ±1 (top and the bottom of the band) leading 
to 



S' '(u>) ~ -T+^r 3 + - + ^r 2 

7T 3 7T 3 
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and 
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Here the first lines in (21,22) correspond to contributions 
from the Fermi points and the second lines to the ones 
from the top and bottom of the band. In particular, 
we see that the only term divergent for u> — > becomes 
exponentially suppressed for temperatures T < J. We 
expect that this picture will remain valid even in the 
interacting case but with possibly renormalized T- and 
w-exponents. This conjecture will be tested numerically 
in the next section. For infinite temperature the integrals 
(19,20) can be expressed as 
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where K, E are the elliptic integrals of the first and sec- 
ond kind, respectively. This confirms the logarithmic fre- 
quency dependence for to — > in this case. 

Our analysis of the interacting case will be based on 
real-time numerical data for Sq\ (t) . In order to calculate 
the Fourier transform we need to extrapolate in time. As 
a guide we will use the long-time asymptotic in the free 
fcrmion case. We can write (18) as 
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For zero and infinite temperature this leads to 29,30 

' \ [J (t) - iH (t)} 2 , T = 

T = oo 



S$'(t) 



\4(t) , 



and 



Sl'{w) 



\[Ji(t)+\H_ x {t)? 



y?(t) 



T = 

T = oo 



(26) 



(27) 



where J„, H n are the n-th order Bessel and Struve func- 
tions, respectively. The long-time asymptotics of these 
functions is given by J n {t) ~ y/2/(irt) cos(t— nn/2 — 7r/4), 
H (t) ~ v /2/(7rt)sin(< - tt/4) + 2/(irt) and #_i(i) ~ 
y/2/(ni) sin(i + 7r/4). It is instructive to derive the long- 
time asymptotics for (25) directly in the T = case. For 
this purpose, consider the integral 
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I(r,t) = - ( dkc l{kr - tcosk ^ 

T J-tt/2 



(28) 



For t > r, the integral has a saddle point at k = 
— arcsin(r/t). That means that I(r, t) for t — > oo is dom- 
inated by contributions from fc ~ 0. These contributions 
can be evaluated by steepest descend methods. For t < r 
and r> 1 the most important contributions come from 
the Fermi points k ~ ±7r/2. Taking both contributions 
into account leads to 



I(r,t) 



— exp[-i(i - tt/4)] cxp[-ir7(2i)] (29) 
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where the first line is the k ~ contribution and the sec- 
ond line the one from k ~ ±7r/2. For the second integral 
in (25) we can do an analogous calculation leading to 
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This result cast some doubt on the field theory results for 
1/Ti discussed in the previous section: The second and 
third term in (30) can also be obtained by bosonization. 
In the interacting case, the exponent of the alternating 
term as well as the amplitudes of both terms will then be- 
come A-dependent. The first term, however, which com- 
pletely dominates for large t, cannot be obtained by these 
methods. In other words, field theory only describes the 
time-dependence of correlation functions for r ^$> 1 and 
t < r. 

By a similar calculation we can also obtain the long- 
time asymptotics at finite temperature, in particular 
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IV. NUMERICS 



Recent progress has made it possible to study the 
real-time dynamics in one-dimensional quantum sys- 
tems by DMRG methods. Here we want to use a 
new variant of the DMRG applied to transfer matrices 
(TMRG) 25 to calculate the longitudinal, real-time auto- 
and nearest-neighbor correlation functions for the XXZ- 
model. The main advantage of this method compared to 
exact diagonalization 31 is that the thermodynamic limit 
can be performed exactly. Therefore time-dependent cor- 
relation functions can be calculated for arbitrary dis- 
tances over a wide temperature range. The time range, 
however, is limited by the fact that the spectrum of the 



reduced density matrix used to truncate the Hilbert space 
becomes dense. In all calculations presented here 400 
states were kept in the real-time TMRG algorithm. A 
detailed analysis of the accuracy of this method and the 
maximum times currently achievable has been presented 
in Ref. 25. In some sense the real-time TMRG method 
is complementary to the calculation of imaginary-time 
correlation functions by standard Quantum Monte-Carlo 
(QMC) 20 ' 21 or TMRG 23 methods: Imaginary-time cor- 
relation functions can be calculated completely due to 
periodicity. The difficulties in this approach arise from 
the analytical continuation which is an ill-posed problem. 
In the real-time domain, on the other hand, the Fourier 
transform is well defined but we have to deal with incom- 
plete numerical data. 

To extrapolate in time we will fit the real part of the 
numerical data for the auto- and nearest-neighbor corre- 
lation function by 

f R (t) =[A + Be-T* cos(fi(t - t ))] /t d (33) 

and the imaginary part by 
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with fit parameters A, B, 7, fl, to, d for the real and 
A, -6,7, ft, Q.2,to,to2,d for the imaginary part, respec- 
tively. These fit functions are motivated by the long-time 
asymptotics in the free fermion case (31,32). Note, that 
7 = 7 = in the free fermion case and that next-leading 
corrections have been taken into account for the imag- 
inary part. The idea is to start for each temperature 
considered with the exactly known fit parameters in the 
free fermion case and then increase the interaction A in 
small steps thus guaranteeing good start values for the 
fit parameters for each anisotropy. As example, the auto- 
correlation function at A = 1 for different temperatures 
and the corresponding fits are shown in Fig. 2. 

We begin our detailed analysis with the case T = 00 
where the correlation functions are real. From the re- 
sults in the previous section for the free fermion point, 
we expect that any singular frequency-dependence will be 
most pronounced in this limit. Particularly interesting in 
this context is how the fit parameter d in (33) evolves as 
a function of anisotropy A (see Table I). The numbers 
obtained clearly show that d au to for the auto- as well 
as djstN for the nearest-neighbor correlation function de- 
crease with increasing A. Furthermore, d au to ~ Ann for 
all anisotropics. This agrees with the expectation that 
the power-law decay for t^>r should not depend on the 
spatial distance r. The autocorrelation function in the 
case T = 00 has been investigated previously on the ba- 
sis of exact diagonalization data for chains up to N = 16 
sites. 31 There, the same fit function (33) has been used 
to analyze the long-time asymptotics and the exponents 
obtained show the same trend as a function of anisotropy 
(see Table I). 
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FIG. 2: Numerical data (solid lines) for the autocorrelation 
correlation function at A = 1 and fits (dashed lines) according 
to Eqs. (33,34). 



TABLE I: The fit parameter d as a function of anisotropy 
A at T — 00 for the auto- and nearest-neighbor correlation 
function, respectively. A variation of the fit region gives an 
error estimate ~ ±0.1 in all cases. For comparison, the values 
dauto obtained in Ref. 31 by a fit to exact diagonalization data 
for a N = 16 site system are shown. 
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The extrapolated numerical data can then be Fourier 
transformed. The results are shown in Fig. 3. At small 
frequencies both correlation functions show a power-law 
divergence Sq\(uj) ~ uj~ a with an exponent a = 1—d and 
d as in Table I. In particular, a ~ 0.3 — 0.4 for A = 1. 
Although this does not agree with the phenomenolog- 
ical theory of spin diffusion by Bloembergen 32 and de 
Gennes 33 which would predict a = 1/2, it is already ex- 
tremely difficult in actual NMR measurements to deter- 
mine whether or not the frequency dependence is singular 
let alone to determine the exponent. From this perspec- 
tive, we might call any kind of divergence a diffusion-like 
behavior. Using this terminology, we conclude that there 
is indeed a diffusion-like contribution to the spin-lattice 
relaxation rate at infinite temperature. Another point 
worth mentioning is the high-frequency tail in Fig. 3 for 
A^0. In the free fermion case all excitations contribut- 
ing to the dynamical structure factor and therefore to 
Sq\(uj) are single particle- hole excitations. The energy 
of these excitations is limited by the bandwidth. For the 
interacting case, however, excitations of multi particle- 
hole type are possible which can carry arbitrarily large 




FIG. 3: So z (w) and S zz (uj) for infinite temperature and dif- 
ferent anisotropics. For A = the exact results (23,24) are 
shown. Insets: Test of the sum rule (35) for each anisotropy. 
The lines denote the exact results, the dots the integrals of 
the numerical data. 



energies. 

A test if the extrapolated real-time data indeed yield 
reasonable results for S§ \ (w) is provided by the sum rule 



Z7T 



(35) 



Because S zz (— ui) — exp(~f3uj)S zz (uj) it is sufficient to 
consider the correlation functions for positive frequencies 
only. For the autocorrelation function S$ z (t = 0) = 0.25 
for all anisotropies and temperatures. For Sl z {ui) and fi- 
nite temperatures, the integrated intensity has to be com- 
pared with numerical data for the static correlation func- 
tion. For infinite temperature, however, Sl z (t = 0) = 0. 
The results of this test are shown in the insets of Fig. 3. 
Next, we consider finite temperatures. Based on the 
exact solution in the free fermion case and on what we 
know from field theory about the contributions from the 
Fermi points, it is reasonable to assume that the ex- 
ponents d in (33,34) are identical and that this expo- 
nent depends on anisotropy only and not on temper- 
ature. We have therefore fixed d — (d au to + rfjvw)/2 



for each anisotropy with d auto , djyw as in Table I. 35 The 
Fourier transform of the extrapolated data for T = 2.0 
and T = 1.0 is shown in Fig. 4 and a check of the sum rule 
(35) in the insets. The results are qualitatively similar 
to the infinite temperature case. In particular, the same 
kind of power-law divergencies for u) — ► and A^ are 
present. Furthermore, a peak in Sq z (u>) starts to develop 
around lo ~ 7r/2. 

For T < 1, however, the divergence at small frequen- 
cies gets strongly suppressed in both correlation functions 
(see Fig. 5) and at T = 0.2 the singular frequency depen- 
dence can no longer be detected. The small kinks visible 
in Fig. 5 for T — 0.2 are not of physical origin. They 
are most likely connected to oscillations in the real-time 
numerical data leading to peaks or dips in Sq\(uj) at the 
corresponding frequencies. Sq z (u>) has been studied pre- 
viously for A = 1 by high-temperature series expansions 
and QMC 21 as well as by a calculation of the imaginary- 
time correlation function using the TMRG algorithm. 23 
The results for this special case (see also Fig. 6(a)) pre- 
sented here are very similar to the ones in these works. 
For example, we also see a peak at w ~ 7r/2 which in- 
creases in height with decreasing temperature and the 
development of a shoulder atw~ 0.7 for T = 0.2. Quan- 
titatively, however, the data in Refs. 21,23 are about a 
factor 2 smaller for all frequencies. As the results here 
fulfill the sum rule (35) for all temperatures with good 
accuracy this suggests that a factor 2 might be missing 
for the numerical data shown in Refs. 21,23. 

The spin-lattice relaxation rates 1/Tf and 1/Tf for 
A = 1 can then be obtained using (11,12) and are shown 
in Fig. 6. The behavior at small frequencies seems to be 
similar in both cases: There is a power law divergence 
oj~ a with a ~ 0.3 — 0.4 at temperatures T > J (remem- 
ber that we set J = 1 here), however, this divergence gets 
strongly suppressed at temperatures T < J. This sug- 
gests that the contributions to the spin-lattice relaxation 
rate with singular frequency dependence behave indeed 
similar to (21,22) found in the free fermion case. 

To analyze the temperature dependence in more de- 
tail, 1/Ti for a fixed frequency ui <C T is shown in Fig. 7. 
The fits in Fig. 7 show that the spin-lattice relaxation 
rate decreases exponentially with temperature in both 
cases and that the scale of the exponential decrease is 
set by J as in the free fermion case. We cannot ana- 
lyze the behavior of 1/Ti at low temperatures in detail 
due to insufficient numerical data. However, the value 
1/Tf (T = 0.2) « 0.31|^| 2 is close to the one predicted by 
the field theory formula (16) if we include logarithmic cor- 
rections to scaling 34 yielding 1/Tf (T = 0.2) w 0.27|A| 2 . 
We also note that 1/Tf only increases by about 30% 
when changing the temperature from T = 0.2 to T = 0.5. 
This is an indication that 1/Tf will indeed be almost 
constant at low temperatures as predicted by field the- 
ory. For the 0(l)-site, our numerical data are consistent 
with 1/Tf -> for T -> 0. 




0.248. 



FIG. 4: Sq z {uj) and Sl z {u) for T = 2.0, 1.0 and different anisotropies. For A = the exact results (19,20) are shown. In the 
insets the integrated intensities (dots) are compared with the exact result S(j z (t = 0) = 0.25 for the autocorrelation and with 
numerical data for the static nearest-neighbor correlation function (lines). 



V. CONCLUSIONS 

The purpose of this article has been to investigate 
if a one-dimensional Heisenberg antiferromagnet has a 
diffusion-like contribution to the spin-lattice relaxation 
rate as has been proposed in Ref. 13 based on 17 NMR 
experiments on Sr2Cu03 . To tackle this problem I found 
it useful to consider the more general XXZ-ca.se which 
interpolates between the exactly solvable free fermion 
and the Heisenberg model we are interested in. For the 
free fermion model I have shown that a contribution to 
the spin-lattice relaxation rate exists which diverges log- 
arithmically for frequency u> — > 0. However, this contri- 
bution comes from the top and the bottom of the band 
and becomes therefore exponentially suppressed at tem- 
peratures T < J. The contributions from the Fermi 
points, on the other hand, do not show any singular 
frequency dependence. I then analyzed the interacting 
case based on real-time numerical data for the auto- and 
nearest-neighbor correlation functions obtained by the 
density-matrix renormalization group applied to transfer 
matrices. 25 The advantage of working in the real-time 



domain compared to imaginary-time methods is that the 
ill-defined analytical continuation of numerical data is 
circumvented. On the flip side, there is no periodicity 
in real-time so the numerical data have to be extrapo- 
lated in time before they can be Fourier transformed. I 
showed that one can do such an extrapolation using the 
long-time asymptotics in the free fermion case as a guide. 
I verified that the results obtained by this method do ful- 
fill the sum rules with good accuracy for all anisotropies 
and temperatures considered. The numerical data for in- 
finite temperature show that the logarithmic divergence 
of Sp* (w) for u> — > in the free fermion case becomes a 
power-law ~ uj~ a in the interacting case. In particular, I 
found a ~ 0.3 — 0.4 for A = 1. With decreasing temper- 
ature these power-law divergencies become exponentially 
suppressed in Sq z (oj) as well as in Sf z (oj). For A = 1, I 
showed that the scale for this exponential suppression is 
still set by J as in the free fermion case. The numerical 
data at low temperatures for 1/Tf (dominated by exci- 
tations with wave-vector q ~ it) as well as the ones for 
1/T-f (dominated by excitations with wave-vector q ~ 0) 
are consistent with the field theory predictions. 
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FIG. 5: So*(w) and Sf z (u) for T = 0.5, 0.2 and different anisotropies. For A = the exact results (19,20) are shown. The 
insets show a test of the sum rule (35). 
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FIG. 6: (a) 1/Tf for the 0(2)-site and (b) 1/7? for the 0(1)- 
site in Sr2Cu03 . 



What does that mean for the NMR experiments on 
Sr2Cu03 ? First, concerning the singular frequency de- 
pendence for u) —> 0, there should be no difference be- 
tween measurements at the copper-, 0(1)- or 0(2)-site. 
Furthermore, a singular frequency dependence should 
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FIG. 7: Temperature dependence of 1/Tf (black dots) and 
1/T? (blue squares) for lu = 0.01 with A = B = 1 
in Eqs. (12,11). The black line is a fit 1/Tf = 0.27 + 
1.62exp(-0.96/T) and the blue line a fit l/T? = 0.0 + 
1.69exp(-0.99/T). 



only show up when the temperature becomes comparable 
to J w 2000 K. A diffusion-like contribution for T « J 
as has been suggested by Thurber et al. 13 based on 17 
NMR measurements at the 0(l)-site cannot be explained 
in a pure Heisenberg model and with the hypcrfine inter- 
action being the only relevant relaxation process. How- 
ever, the evidence presented in favor of such a contribu- 
tion is rather weak. There is no reason to assume that the 
spin-lattice relaxation rate in the limit of infinite mag- 
netic field h is given by the field theory result where the 
effect of the magnetic field on the spin-spin correlations 
has been ignored. In fact, the magnetic field can only be 
ignored if T 3> \h\. Without this limiting value, however, 
the data in Fig. 3(d) of Ref. 13 are also consistent with 
having no frequency dependence at all. If that is the 
case, the only part of the data which is not consistent 
with the simple Luttinger model picture is the tempera- 
ture dependence of the relaxation rate at the 0(l)-site, 
1/{T?T) ~ const + T compared to 1/(T?T) ~ const ex- 
pected from field theory. As contributions from q ~ ir 
are strongly suppressed, this next-leading temperature 
dependent term most likely has to do with corrections to 
the simple <5-peak (14) obtained from the Luttinger model 
for the dynamical structure factor at small q. From re- 



cent studies at zero temperature we know that a finite 
band curvature will broaden the <5-peak and lead to in- 
teresting singularities at the lower and upper thresholds 
as well as to a high-frequency tail. 8-10 At finite temper- 
atures, spectral weight will also appear below the lower 
threshold. Further research, if these corrections can in- 
deed explain the measured temperature-dependence of 
l/(T^T) is necessary. Finally, I want to remark that in 
the free fermion case the next-leading term at low tem- 
peratures is of order T 3 (see Eqs. (21,22)) and not T 2 . 
That suggests that if a T 2 -contribution exists, it should 
either have an amplitude which vanishes for A = or 
otherwise the exponent has to change as a function of A. 
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